SYMMETRICAL AIRFOIL GENERATION PROGRAM, JEF RASKIN, 1984 


While this program is nothing that experienced programmers will learn much 
about coding from, it demonstrates the power of being able to embed executable 
code into the text produced by a word processor. What this does is to make 
the mechanics of documentation simple. Using a typical program editor for 
extended text is usually bothersome enough a chore so that our programm comments 
become terse and more difficult to read. I have observed that when I document 
a program in this thorough fashion, it often runs on the second or third try, 
whereas when I do not document carefully, much more debugging is required. 


please note that this text you are reading is part of the program (or, 
equivalently, that the program is part of the text). 


The program itself plots, on an Apple II, given Information Appliance's 
Swyftware environment, a family of airfoils that I am developing for my 
aerobatic model airplanes. It does it by a transformation of a circle, such as 
the Joukowsky transformation does. This permits relatively easy calculation of 
the lifting characteristics of the airfoil, although a program to do that is 
not shown here. 


The program begins by setting up graphics mode and a few constants. 
First, for graphics mode: 


90 HGR 
The constants include 
100 pi = 3,.1415926 


and one that determines how much the airfoil will "cusp" or be concave toward 
the rear. 


110 c = 0,81 
Another important variable controls the thicknes of the airfoil. 
120 £ = 0.22 
The program works by using the usual pair of formulaes 
x =r cos t y=rsint 
to generate a circle parametrically by letting it vary from 0 to 2 * pi in 


convenient steps. r is chosen to be 0.5 so that a circle of diameter 1 results. 


t= 0 
0.5 * cos (t) 
0.5 * 


The circle is centered on the origin, and we first move it to the right by 0.5 
so that all the x-values are positive, and are, in fact, in the range (0,1). 


160 x = x + 0.5 


To understand how this transforms a circle into an airfoil rounded at one 
end and pointed at the other, note that if we multiply each y-value by the 
corresponding x-value we will, near x=0, have very small y-values, yielding a 
point. But as we approach x=1 the curve approaches a circle, 


200 y = y * x*c 


The exponent of x changes x's effect. Many functions of x will do for 
this role. 


Airfoils tend to be long and thin. Thus we want to squash the shape in 
the y-direction. The thickness ratio is thus applied to y. This could have 
been done in the previous line of code, but it is separated out for clarity. 


aloy=yk#et 


We now have -f < y < f and 0 <= x <=1. This has to be translated to 
Apple screen coordinates where 0 <= x < 280. Again, these transformations 
could have been included in some of the earlier program statements, but for 
clarity the plotting is separated from the generation of the foils. The 250 is 
chosen as conveniently large for the screen. This will make x vary from 0 to 
250 


400 x = x * 250 


and to keep things square, we apply the same scaling to y, remembering that y 
has already been decreased by the f factor. We also have to make y always 
positive. Obviously, we will go off-screen if f is too large, but we will 
assume smart users who read this documentation and behave accordingly rather 
than bother with too much error checking. 


410 y = y * 250 + 100 


Now to plot these points, and close the loop, and then go into an infinite 
loop so that the computer will not mess up the screen with messages. 


500 HPLOT x, 
510 NEXT t 
520 GOTO 520 


All that's left is to type and use the CALC command on ‘runt. You will 
have to modify the program ta run on systems other than SwyftWare-equipped 
Apple II's. 


